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ABSTRACT 

Synchrotron radiation of ultra-relativistic particles accelerated in a pulsar wind neb- 
ula may dominate its spectrum up to 7-ray energies. Because of the short cooling 
time of the 7-ray emitting e*, the 7-ray emission zone is in the immediate vicinity of 
the acceleration site. The particle acceleration likely occurs at the termination shock 
of the relativistic striped wind, where multiple forced magnetic field reconnections 
provide strong magnetic fluctuations facilitating Fermi acceleration processes. The ac- 
celeration mechanisms imply the presence of stochastic magnetic fields in the particle 
acceleration region, which cause stochastic variability of the synchrotron emission. 
This variability is particularly strong in the steep 7-ray tail of the spectrum, where 
modest fluctuations of the magnetic field lead to strong flares of spectral flux. In par- 
ticular, stochastic variations of magnetic field, which may lead to quasi-cyclic 7-ray 
flares, can be produced by the relativistic cyclotron ion instability at the termination 
shock. Our model calculations of the spectral and temporal evolution of synchrotron 
emission in the spectral cut-off regime demonstrate that the intermittent magnetic 
field concentrations dominate the 7-ray emission from highest energy electrons and 
provide fast, strong variability even for a quasi-steady distribution of radiating parti- 
cles. The simulated light curves and spectra can explain the very strong 7-ray flares 
observed in the Crab nebula and the lack of strong variations at other wavelengths. 
The model predicts high polarization in the flare phase, which can be tested with 
future polarimetry observations. 

Key words: shock waves — turbulence — ISM: supernova remnants — gamma-rays — 
supernovae: individual (Crab nebula) 



1 INTRODUCTION 

Strong flares of a few days duration have been discov- 
ered recently by the AGILE and Fermi 7-ray observatories 
in the Crab nebu l a at energies above 100 MeV (see e.g . 
iTavani et all I2OIII : lAbdo. et af] I2OIII : IVittorini etaj l201ll . 
and references therein). The most striking features of the 
flares are the extreme amplitude of the photon flux changes, 
especially at energies above the exponential cut-off energy 
of the quiescent spectrum, and the fast hour-timescale vari- 
ability. While the exponential cut-off energy Ec of a typical 
quiescent 7-ray spectrum is ~ 100 MeV, the cut-off energy 
of > 500 MeV was found in the April 2011 flare spectrum. 
This value of Ec exceeds the energy Ec ~ nieC'^ /a ~ 100 
MeV (where a — e^/hc is the flne-structure constant), con- 
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sidered to be the maximal cut-off energy in the synchrotron 
models of 7-ray emissi on from the Crab and other pulsar 
wind nebulae (PWNe) Jde Jager et al]| 19961 : lUzdenskv et all 
I2OIII : IStriani et al.ll201ll '). The value Er is the synchrotron 
photon energy emitted by an electron whose energy is such 
that the synchrotron cooUng time is equal to the charac- 



teristic gyration time uig ^ (see e.g. Guilbert et al.l Il983l : 



Ide Jager et~al]|l996l : lAtovan fc Aharonianlll996l '). The par- 
ticle gyration time is considered to be the fastest accelera- 
tion time in a plasma system with frozen-in magnetic fleld 
of the r.m.s. amplitude {_B^)^/^ that exceeds the electric 
fleld magnitude. The value of Ec corresponds to the elec- 
tron Lorentz factor 7™ that can be found from the equa- 
tion £syn{'~tm) = face (7m). Both the syuchrotrou loss rate 
£syni'y) and the electron acceleration rate face (7) depend 
on the moments (often just (-B^)) of the stochastic magnetic 
fleld. 
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Figure 1. Light curves of synchrotron emission at 5 keV (dot- 
dash Une), 500 MeV (dashed hne), 1 GeV (dotted hne) and 2 
GeV (short-dash line) as a response to an imposed fluctuation 
with magnetic field B{t) (solid line) simulated in Model I. The 
light curves are normalized to maximal intensity. The background 
magnetic field in the emission region was modeled as a stochastic 
gaussian field of (B2(0))1/2 = 0.2 mG. The imposed fluctuation 
B{t) (solid line) is localized in a stripe of a O.OIAq thickness 
(> rf). The maximum of B{t) = 1 mG is at ci/Ao = 0.6. The 
spatial scale Aq ~ 2 x 10^® cm for the photon energy E ^ 1 GeV 
provides the photon variability timescalc of about 10"'' s in the 
GeV regime. 



In this work we consider the case when the formation 
length, rf = nie /e (B^)^^^ , of incoherent synchrotron ra- 
diation is much smaller than the typical synchrotron cool- 
ing and acceleration lengths, while the typical wavelength 
A of the fluctuating magnetic field is larger than r/. In this 
case Ssyni'ym), £acc{'~fm), and Ec are determined by the same 
r.m.s. value of the fluctuating magnetic field, and the bulk of 
the relativistic electron distribution may vary on scales that 
are much larger than the gyration radius rg — ^rf. Since 
the synchrotron emissivity of a power-law electron distribu- 
tion with spectral index p is proportional to 5(^+1)/^^ the 
local emissivity sharply grows with B for large p values. This 
means that the synchrotron radiation in the cut-off regime 
(which corresponds to large effective p values) is governed 
by high statistical moments of the stochastic magnetic field 
distribution, and it is intermittent. The intermittency ef- 
fect implies that rare strong peaks of the magn etic field dis- 
tribution dominate the synchrotron emission (| Bykov et al.l 
[2OO8, 2009). It is particularly important in the synchrotron 
cut-off regime, when the typical size of the distribution of 
radiating electrons (the synchrotron cooling length) can be 
comparable with the correlation length of strong magnetic 
field fluctuations. For instance, this is expected to be the 
case in supernova shells, where magnetic fluctuations are 
produced by instabilities of anisotropic distributions at the 
max imal energy of par ticles accelerated in the source (see, 
e.g.. lBvkov et al.ll201ll . and references therein). 

Since the source emission in the cut-off regime is domi- 



Figure 2. Normalized spectra of synchrotron radiation at two 
different time moments ct/Ag = 0.2 (solid line) and 0.6 (dashed 
line), which model the quiescent and flare spectra, respectively 
(see Fig.[T]l. The dotted curve shows the contribution of the vari- 
able magnetic field. The power emitted in the GeV fiare is about 
2x10'^'' erg s~^, for the Crab parameters. 



nated by just a single (or a few strongest) concentration(s) 
of the stochastic magnetic held, the light curve of the source 
in this regime reflects the lifetime of the magnetic concentra- 
tions rather than the electron acceleration/losses timescales. 
Fast temporal variations will appear even for a quasi-steady 
electron distribution. The light curve and the spectral be- 
haviour in the synchrotron cut-off regime are determined by 
statistical characteristics of the magnetic field, which can be 
described by the probability distribution function (PDF) of 
magnetic fluctuations, P{B). 

It is worthwhile to note that intermittent magnetic 
flelds can be found in quite different circumstances. For in- 
stance, non-Gaussian distributions of fluctuations, exhibit- 
ing gradual tails at large field amplitudes, have been found 
in the Earth magnetotail after the current disruption asso- 
ciated with magnetospheric substorms observed by Geotail 
and Cluster satellites. The energy injection during the sub- 
storms feeds an energy cascade to small-scale fluctuations 
with the correspondi ng increase of intermittency (see, e.g., 
IZimbardo et al.(l2010l . for a recent review of magnetospheric 
observations). 

In this Letter we show that the model of synchrotron 
emission in fluctuating magnetic fields, with account for the 
intermittency in the spectral cut-off regime, can explain the 
nature of the 7-ray flares observed in the Crab. We demon- 
strate the effect of intermittency on the GeV regime emission 
using two models. In first model, we simulate the spectra of 
accelerated electrons and positrons in a simple kinetic model 
of diffusive Fermi acceleration in the termination shock of 
striped wind with account for synchrotron losses. Then we 
construct the GeV regime flare light curve and spectra by in- 
tegrating the synchrotron emissivity of spatially inhomoge- 
neous particle distribution in the shock downstream with im- 
posed magnetic held variation. Second model demonstrates 
the effect of the magnetic held PDF shape on the syn- 
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chrotron photon spectra in the cut-off regime. The emission 
is produced by Fermi-acceierated pairs in the spectral cut-off 
regime, in which the acceleration is balanced by synchrotron 
losses. 



2 MODELING 

The synchrotron origin of the observed gamma-ray flares 
assumes the presence of pairs accelerated to PeV energies. 
The thickness of the 7-radiating region depends on the cool- 
ing rate in the magnetic field: A ~ TsynC ~ 1.5 x 10^^ ■ 



B. 



-3/2 



-1/2 



cm, where Bmc is the r.m.s. magnetic field 



in mG, and EgcV is the photon energy in GeV. The size 
of the 7-ray emitting region is very small compared to the 
size of the nebula. The thickness of the layer A is estimated 
using the standard electro n synchrotron cooling time (e.g. 
iRvbicki fc Lightmanlll97^ ) of the electron radiating photons 
at the peak of the power spectrum of synchrotron emis- 
sion. The particle acceleration mechanism in PWNe is not 
fully understood yet, although some bas ic features and con- 
straints have been e stablished (see, e. g., 'Kenn el fc Coronitil 
1984 lArond I^^P^' I2OIII: iKeshet et al. 2009j; iKirk et all 
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Bvkov fc TreumannlboTiriUzdenskv et aLlbOlll v 

Recent models of particle acceleration in PWNe con- 
sider a relativistic wind in the equatorial plane with toroidal 
stripes of o pposite magnetic field polar ity, separated by cur- 
rent sheets. ISironi fc Spitkovskvl (|201ll ) modeled particle ac- 
celeration and magnetic field dissipation at the termina- 
tion shock of a relativistic striped wind using 2D and 3D 
particle-in-cell (PIC) simulations. They found a complex 
structure of the flow in the vicinity of the termination shock. 
In that model, the shock-driven reconnection in the down- 
stream transfers the magnetic energy of alternating fields of 
the striped wind to accelerated relativistic pairs. The en- 
ergy spectra of electrons accelerated in the reconnection re- 
gion take the form of power laws, A''(7) (x 7"*^, with spec- 
tral indices p ~ 1.5 that ma tch the radio-optical observa- 
tions of the Crab nebula fe.g. iBietenholz et al.|[l997l : iHesteJ 
l2008l : lArendt et aill201ll 'l. The accelerated particles can es- 
cape ahead of the shock and generate magnetic fiuctuations 
in the upstream by the filamentation and/or Weibel type 
instabilities. The turbulence generated by the instabilities 
can alleviate Fermi-like diffusive process that accelerates X- 
and 7-ray emitting electrons. This is an important finding 
as it simultaneously addresses two problems - the termina- 
tion shock formation in magnetized PWN winds and particle 
acceleratiorQ. Therefore, the combined action of the recon- 
nection processes and shock acceleration is expected in the 
intense equatorial pulsar wind. Extensive PIC simulations 
in a wide dynamical range are needed to demonstrate the 
feasibility of this approach and identify the mechanism of 
particle acceleration to PeV energies and their synchrotron 
emission. However, to construct the synthetic spectra of the 



^ Note that the standard model of diffusive Fermi accelera- 
tion in a transverse relativistic shock in a non-striped uniform 
wind encounters problems when appli ed to PWN terminatio n 
shocks, see, e.g.. iNiemiec et al.l 1I2OO6I I; iPelletier et"ai] | |2009|) : 
iBvkov &: TreumannI 1I2OIII') . 



Crab nebula at 7-ray energies (where the synchrotron cool- 
ing is very important), one needs to account for the radiative 
reactio n force, which is not y et attainable i n the PIC sim u- 
lations (|Sironi fc Spitkovskvl [2009: Nishik awa. et al.ll201ll '). 

Kinetic models can be used to simulate particle ac- 
celeration due to the repetitive interaction of electrons 
with magnetic turbulence in the energetic outflows with 
acc ount for the synchrotron cooling. The kinetic model 
by iBvkov fc Meszaro^ (| 19961 ) for particle acceleration by 
both relativistic and transrelativistic shocks, accompanied 
by broad dynamic spectra of magnetic fluctuations with vi- 
olent motions of relativistic plasma, predicts a hard broken 
power-law electron distributions with slopes 1 ^ p ^ 2. In 
this model, particles are accelerated by strong magnetic fluc- 
tuations on timescales comparable to their gyration period 
in the r.m.s. magnetic field. 

In the simulations described below we complement the 
kinetic model for particle acceleration in the vicinity of the 
striped wind termination shock with the synchrotron cooling 
effects to account for the spectral cut-off regime. To estimate 
the spectra of nonthermal leptons accelerated downstream 
of the PWN termination shock by Fermi mechanism, one 
can use a Fokker-Planck-type kinetic equation: 



dN 1 a 2 

oz 7^ 



97 



dN 

D{f) — +a{'y)N 



(1) 



where 2 is the coordinate along the shock normafl This 
equation, averaged over the ensemble of strong electromag- 
netic fluctuations in the vicinity of the wind termination 
shock, accounts for diffusion and advection of electrons in 
phase space due to interactions with the fluctuations. The 
term with the momentum diffusion coefficient D{'y) corre- 
sponds to the stochastic Fermi acceleration, k{'y) is the fast 
particle spatial diffusion coefficient, u is the fiow velocity 
component along the shock normal, and 0(7) is the energy 
loss rate of an electron due to synchrotron radiation. 

Model I. To construct the synchrotron emission spectra 
and flaring light curves in the diffusive shock acceleration 
model, we simulated spatially inhomogeneous accelerated 
pair distribution downstream of the termination shock of 
the striped wind using Eq. [T] Short-scale magnetic fluctu- 
ations are required to be present upstream of the shock to 
allow an efficient diffusive Fermi acceleration in the trans- 



verse relativistic shocks (s ee e.g. IBvkov fc Treumanr]|20l"ll . 
and footnote ^). Recentlv. ISironi fc Spitkovskv 1 20111 ) have 



found in the PIC simulations that the fiuctuations can be 
generated upstream of the striped wind termination shock. 
We assume the fluctuations provide the Bohm diffusion with 
k{"/) ~ erg (7)73. The stochastic Fermi acceleration was ne- 
glected in the model (i.e., D = 0). 

To illustrate the intermittency effect in the cut-off 
regime, we simulated a light curve and spectra for a mag- 
netic field fiuctuation imposed in the GeV photon emitting 
region of scale size Ao — 2 x 10^^ cm (for a quiescent 
magnetic field Bo = 0.2 mC, downstream of the termina- 
tion shock). The fluctuation 5B{t) is localized in a stripe 
of a O.OlAo width and has the time dependence shown by 

^ Since the scale size A of the PeV electron distribution is much 
smaller than the termination shock radius, the problem can be 
considered as one-dimensional. 
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the solid line in Figure [T] Such magnetic field variations 
may be produced by the relativistic ion cyclotron instabil- 
ity at the termina tion shock. The instability, proposed by 
ISpitkovskv fc Arons ( 200 4) to explain the origin of the opti- 
cal wisps in the Crab nebula, can produce quasi-cyclic 7-ray 
flares in our model. The scale Ao, used in our simulations, 
correspon ds to about 0.1 of t h e mag netic field limit cycle 
found by ISpitkovskv fc Arond (|2004l i (see their Figure 2). 
The light curves of the 7-ray (0.5 , 1 GeV and 2 GeV) and 
X-ray (5 keV) emission show the strong response of the 7- 
ray emission in the cut-off regime, while the response is very 
modest at X-ray energies. Note that the variability time in 
the GeV regime is shorter than the imposed field fluctuation 
and the cooling time of PeV electron distribution. This is be- 
cause the emission in the cut-off regime is governed by high- 
order momenta of the magnetic field. The energy loss rate 
0(7) that determines the electron cooling depends on (B^). 
Therefore, the stochastic magnetic field realizations in the 
emission region with the same {B^) but different high-order 
momenta (determined by their PDFs) would correspond to 
the same particle distributions. However, their photon spec- 
tra in the cut-off regime are very different. This effect is 
clearly seen in Figure [3] discussed below. The synchrotron 
flares in the cut-off regime will appear even in the case of 
steady electron distributions. 

In Figure [2] we show the simulated spectra of syn- 
chrotron emission, generated downstream of the termina- 
tion shock, for 77 ~ 0.5, corresponding to quiescent regime 
(solid line) and fluctuation maximum (dashed line). They 
are similar to the quiescent and flare spectr a observed in 
the Crab nebula by AGILE and Fermi (se e iTavani et al.l 
I2OIII : lAbdo. et al]|201ll : IVittorini et"aill201ll . and references 
therein) . 

Model II. Apart from the diffusive shock acceleration 
a variety of other particle acceleration mechanisms can be 
important in the region. The coalescence of magnetic is- 
lands, particle reflection by magnetic islands (both first 
and second order Fermi- type processes), are expected to 
be in action, as it occurs in the Earth magnetosphcrc (e.g. 
Drake ct al. 2006 ; Zolcnyi ct al. 2010; Siroru & Spitkovsky 
20 111 : iUzdenskv et al.l I2OIII : iDaughton et all I2OIII ). When 



the source of strong turbulence is quasi-steady on timescales 
longer than aj~^(7), a simple analytical treatment of the 
problem is also possible. For the case of fast stochastic Fermi 
acceleration (comparable to the particle gyration period) by 
an energetic plasma outflow with strong magnetic turbu- 
lence, the momentum diffusion coefficient in Eq[T] takes the 
form D{pj) — 7^ rju>g{^), where 77 < 1 is the particle acceler- 
ation time measured in the gyration times. Note that both 
ujgi-y) = e{B2)i/2/(^^c7) and 0(7) = Arl{B^)-i^ / [Qniec) 



(where ro 



/me 



depend on the same ensemble- 



averag ed value (B^). 

The asymptotical shape of the particle spectrum in the 
cut-off regime is 



A''(7) oc 7 exp 



d7a(7)/D(7) 



(2) 



It is important that, while the index p depends on the 
turbulence spectrum and system geometry, the exponen- 
tial cut-off in the particle spectrum is rather universal: 
7V(7) oc 7-f exph(7/7o)'], where -yl = 9er?/(2rg {B^)^'^), 
i.e., 70 ~ 5 X 10^ »7 (-B^g)"^''*. The synchrotron emissivity 



e[ui, B) in a local magnetic field B, which is assumed to be 
uniform on spatial scales larger than r/ — m^c^ /eB, is given 
by the equation 



e{uj,B,z) 



27rmc2 



d'y'y^N{z,'y)R{iu/Lj,), 



(3) 



where tOc — 3eB'y^ /2meC is the characteristic frequency of 
synchrotron radiation. Approximate analytic expressions for 
the fu nctio n R(x) were derived by Crusius fc Schlickeiserl 
(|l986h and lZirakashvili fc AharonianI l|20o'7[ r The spectrum 
of synchrotron emission from the downstream region filled 
with strong magnetic field fluctuations can be expressed as 



dB dze{oj,B,z)P{B). 



(4) 



To illustrate the effect of the magnetic field fluctuations 
in the cut-off regime, where high-order statistical moments 
dominate the integral in Equation (j4}, we used the PDF of 
magnetosonic type fluctuations, which corresponds to the 
wisp stru ctures s e en in the polarized optical images pre- 
sented bv iHesterl (|2008^. The PDF has the form P{B) = 
C„ exp(-6"/e„), with n = 1 and 2, and b = \B - Bo\/Bo. 
Here Cn is the normalization constant, and 0„ is the di- 
mensionless width of the distribution. The simulated syn- 
chrotron spectra are presented in Figure [3] for the gaussian 
(n — 2), exponential (n = 1), and non-fluctuating magnetic 
fleld distributions. The characteristic frequency ojq for the 
Lorentz factor 70 is given by hujo ~ 277? m^c'^/ (4q:) ~ 470 77 
MeV. The synchrotron curves in Figure [3] are simulated for 
6n = 1 (i.e., for the case of strong fluctuations). The results 
illustrate a strong effect of the PDF shape on the spectral 
behaviour in the cut-off regime (even for a fixed rms field 
{B'^)), in contrast to a very modest effect in the power-law 
regime at a; < uiq. Thus, Model II demonstrates that a re- 
construction of the PDF tail of magnetic fluctuations in the 
synchrotron emission region (with the same (B^)) would re- 
sult in a strong change of the synchrotron photon spectrum, 
similar to that observed in the GeV flares in the Crab neb- 
ula. 



3 DISCUSSION 

In Model I, 7-ray photons are radiated downstream of the 
PWN termination shock by ultrarelativistic electrons. The 
particles are accelerated by some kind of Fermi mechanism 
in the vicinity of the termination shock in the striped wind 
with reconnecting magnetic flelds. Electrons in the sharply 
decreasing high-energy tail of the distribution function radi- 
ate synchrotron 7-ray photons in fluctuating magnetic flelds, 
which results in flaring behaviour in the cut-off regime (i.e., 
at the high-energy end of the synchrotron spectrum). 

The compressed magnetic fleld structures, moving 
through the very narrow 7-ray emitting layer just down- 
stream of the termination shock, result in 7-ray flares 
with the strong spectral variations shown in Figure [2l The 
timescale of the 7-ray photon variability in the flare is a 
fraction of the time Ao/c ~ 5 x lO'' s in the cut-off regime 
because of the effect of high-order statistical moments. At 
the same time, the variability at energies below the cut- 
off energies is slower, and of much smaller amplitudes, as 
it is clearly seen from comparison of the light curves at 
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Figure 3. Surface spectral density of synchrotron radiation 
in the cut-off regime (in arbitrary units) simulated for different 
probability distribution functions of fluctuating magnetic field in 
the source (Model II). The characteristic frequency is defined by 
huo Si 470 ri MeV. 



5 keV and above 500 MeV in Figure [T] This is consistent 
with the lack of strong variatio ns in t he Crab at the X-ray 
and optical energie^ (jTennant et all bOllI : [C araveo et al. I 
The frequency of the 7-ray flares is determined in 
our model by the frequency of occurrence of extreme values 
of magnetic fluctuations. The strong magnetic field compres- 
sions could be associated with the variable wisps o bserved 
with t he quasi-periodic cycle of a year timescale (e.g. lScargld 
1 19691 : iHester et all I2OO2I : iBietenholz et al.ll2004). To model 
the dynamics of the wisps, ISpitkovskv fc Arond l(2004l ) simu- 
lated outward-propagating magnetosonic waves downstream 
of the shock generated by relativistic cyclotron instabil- 
ity of gyrating ions. The process achieves a limit cycle, in 
which the waves are launched with periodicity on the or- 
der of the ion Larmor time (a year timescale). Compres- 
sions in the magnetic field and pair density, associated with 
these waves, can reproduce the behaviour of the wisp and 
ring f eatures. The high-re solution relativistic MHD mod- 
els bv lCamus et al.l (|2009l ) have revealed a highly variable 
structure of the pulsar wind termination shock with quasi- 
periodic behaviour within the periods of about 2 years, and 
MHD turbulence on scales shorter than 1 year. 

The synchrotron emission in the 7-ray regime is domi- 
nated by infrequent quasi-coherent structures, and therefore 
it should be highly polarized, that can be tested with future 
time resolved polarimetry observations. 

The strong 7-ray variability predicted by our model is 
consistent with that observed in the Crab nebula. No simul- 
taneous strong flares at other wavelengths is expected in our 
model. The 7-ray flares can occur even for a steady electron 



^ Note that we are discussing variations in the flux integrated 
over the downstream region (over 2 in our ID model). The lack 
of strong variations in the flux does not contradict local variations 
in the synchrotron intensity, such as the moving wisps. 



distribution function with the maximal energies not exceed- 
ing the limit of particle acceleration by the Fermi mechanism 
with strong magnetic fluctuations at the extended wind ter- 
mi nation shock . This is the main difference from the model 
bv lYuan et al.l (|201ll ). which explains the flares by varying 
maximal electron energies in isolated knots that must trans- 
fer a substantial power (above 10^® erg s~^) to the observed 
GeV photons. The variability of the maximal energies of the 
electrons would result in simultaneous GeV and 100 TeV 
regime flares since the electrons producing the synchrotron 
GeV photons radiate also 1 00 TeV regime photons b y inverse 
Compton scattering. Also, iBednarek fc Ide3 (|201ll ) pointed 
out that the TeV photons produced by inverse Compton 
scattering of soft radiation by the variable distribution of 
accelerated electrons should be variable on timescales simi- 
lar to those observed at GeV energies by AGILE and Fermi- 
10" LAT. 

Recently, iKomissarov fc Lvutikovl (|201ll ) proposed a 
model in which the Crab flares originate in the inner knot 
(within about 1" from the Crab pulsar) with strong Doppler 
beaming effects. In our model the region of GeV photon 
emission is about 10" away from the pulsar. The Doppler- 
boosted s ynchrotron emission fr om a corrugated shock, pro- 
posed bv iLvutikov et all l|201ll ). would be accompanied by 
flares in the cut-off regime of the inverse Compton TeV pho- 
tons. This is different from the prediction of our model, in 
which the strong variability of GeV synchrotron emission is 
due to magnetic field variability, strongly amplified in the 
cut-off regime, while the amplitude of TeV emission varia- 
tion is expected to be much less prominent. 
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